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LARGE SPACE STRUCTURE CONTROL PROBLEM 


Large Space Structures (LSS) exhibit characteristics which make the LSS 
control problem different from other control problems. LSS will most likely 
exhibit low-frequency, densely spaced and lightly damped modes. In theory the 
number of these modes is infinite. Because these structures are flexible, Vi- 
bration Suppression is an important aspect of LSS operation. There are a num- 
ber of implementability issues which must be dealt with by any "space realiza- 
ble" actuation and sensing scheme. In terms of Vibration Suppression, we would 
like the control actuators to be as low mass as possible, have infinite band- 
width, and be electrically powered. In addition, we argue that actuators which 
produce "internal forces" in the structure have distinct advantages for the 
Vibration Suppression application. Since velocity sensing may be very dif- 
ficult at low vibration levels and low frequencies, we prefer to use strain 
as the only measurement. Finally, we propose that actuators be built into the 
structure as dual-purpose structural elements in the interest of efficiency of 
design. 


• LOW FREQUENCY, DENSELY SPACED AND LIGHTLY DAMPED MODES ARE COMMON. 
ACCURATE SYSTEM CHARACTERISTICS ARE DIFFICULT TO OBTAIN 

• VIBRATION SUPPRESSION IS AN IMPORTANT ASPECT OF LSS OPERATION 

• ACTUATOR REQUIREMENTS ARE LIGHTWEIGHT, INTERNAL FORCE PRODUCING 
ELECTRICAL POWERED, INFINITE BANDWIDTH, etc. 

• ACCURATE VELOCITY SENSING MAY BE UNREALISTIC 

• STRAIN SENSORS SHOULD BE CONSIDERED 

• INTEGRATED DUAL PURPOSE LOAD CARRYING/ACTUATION MEMBERS SHOULD 
BE CONSIDERED 
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STIFFNESS CONTROL 

Initial work investigating vibration suppression in LSS using internal 
forces centered on the one-dimensional vibrating string. The string has low 
inherent out-of -plane stiffness, like some Large Space Structures. It was 
found that by varying the tension in the string as a function of state vari- 
ables and time, damping could be introduced. The lower figure shows a plot of 
the motion of the string in the phase plane. The distance of the curve from the 
origin is an indication of the energy in the motion at a particular point in 
time. The damping is evident by the spiraling of the locus into the origin. 

VIBRATING STRING 



d 2 y = _3 

P dt 2 


T = T y 


= J-/ T iv\ 

3x \ 3x ) 


e = +0.2 


0.028 LINEARIZED SOLUTION 



EXAMPLES OF ACTIVE MEMBERS 


The concept of an active member is to replace a passive structural ele- 
ment , such as a diagonal of a space-truss beam, with a structure which is also 
a control actuator and sensor. We propose a piezoelectric active member for 
the control of LSS. Such devices would consist of a piezoelectric actuator 
and sensor for measuring strain, and screwjack actuator in series for use in 
quasi-static shape control. Several concepts for active-members are shown. 
One variation is to beam a laser through a hollow strut to measure movement 
between the two ends of the member. We envision these devices as being self- 
contained, possibly containing their own electronics for effecting Vibration 
Suppression. 
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PIEZOELECTRIC 

SENSOR/ACTUATOR 


SCREW JACK 


A = ACTUATOR 
S = SENSOR 




TRUSS MODEL 


ACTIVE MEMBER WITH LASER 
DISPLACEMENT SENSOR 




FEASIBILITY STUDY — PIEZOBEAM EXPERIMENT 


In order to investigate the feasibility of using piezoelectric active mem- 
bers to perform Vibration Suppression in LSS, a simple experiment was design- 
ed. The objective of the experiments is to simulate an active member using 
piezoelectric ceramic thin sheet material on a thin, uniform cantilever beam. 
The structure was designed to have low stiffness, low mass density, and to have 
a first mode at 5 Hz . We use collocated piezoelectric ceramics as both actua- 
tors and strain sensors. The layup of the ceramics and the dimensions of the 
composite piezobeam are shown. 


• OBJECTIVE - SIMPLE DEMONSTRATION OF VIBRATION SUPPRESSION 

• INTEGRATED STRUCTURAL MEMBER/ACTUATOR 

• LSS CHARACTERISTICS 

• LOW STIFFNESS 

• LOW MASS 

• HIGH MODAL DENSITY AT HIGHER FREQUENCY 

• COLLOCATED ACTUATORS/SENSORS 

• SPACE REALIZABLE APPROACH 



LENGTH: 

WIDTH: 

THICKNESS: 

MATERIAL: 


12.50 in 
0.648 in 
0.020 in 
ALUMINUM 


1.25 in 1.25 in 

0.50 in 0.25 in 

0.0095 in 0.0095 in 

LEAD-ZIRCONATE-TITANATE <PZT) 
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TEST SET-UP 


The cantilever beam was supported in a vibration test fixture shown in the 
figure. The beam was supported in a clamping flange which was bolted to a lin- 
ear bearing table. The table was excited by means of a stinger attached to a 
small shaker. A wide variety of waveforms were used to test the open-loop and 
closed-loop performance of the piezobeam. 
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MECHANICS OF PIEZOELECTRICS 


The piezoelectric ceramic material is an inherent electromechanical trans- 
ducer. If an electric field is applied to the material, it tends to strain by 
an amount proportional to the strength of the applied field. The proportion- 
ality constant is the d 31 coefficient. If, on the other hand, the material is 
stressed, an electric field is generated spontaneously. The proportionality 
constant between stress and generated electric field is the g 31 coefficient. 
Both the d 31 and the g 31 coefficients are material properties of the piezoelec- 
tric . 

The piezoelectrics are arranged on the test beam in a sandwich fashion. 

The actuators are arranged such that a voltage applied to the outer electrode 
surfaces causes one ceramic to expand while the other contracts . Since the ce- 
ramics are adhered to the beam, a bending moment is produced. Similarly, the 
bending of the beam stresses the sensor ceramics which in turn produce a volt- 
age which is measured. 


(a.) ACTUATOR PIEZOELECTRIC 

Ef POLARITY 

ceramic yAzzzzzazzzzzz. z 
h 10 H 

(b.) SENSOR PIEZOELECTRIC 

E f POLARITY 

° -*v fyzzzz%zzzzza - ° 


£ E“ d 31 E f< 


Ef - * 931 °- 


V a = ACTUATOR VOLTAGE 

d 31 = PIEZOELECTRIC STRAIN CONSTANT 

g 31 = PIEZOELECTRIC VOLTAGE CONSTANT 

t a = THICKNESS OF THE ACTUATOR CERAMICS 


PUSH-PULL BENDING ACTUATOR 


CERAMIC i 
CERAMIC 



if 

D 


ALUMINUM 
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POLARITY 
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MECHANICS OF PIEZOELECTRICS (CONT.) 


The moment applied by the piezoelectrics is determined by integrating the 
stress produced. The magnitude of the applied moment is found to be propor- 
tional to the width of the actuator ceramic W a , the product of the Young ’ s Mod 
ulus and piezoelectric strain constant d 31 E a , the "lever arm" (or distance 
from the neutral axis) , and the applied voltage V a . The measured voltage was 
about 25°/, less than the predicted value which is consistent with the simplify- 
ing assumptions of the analysis . 


ACTUATOR INDUCED BENDING MOMENT 



M = f o x (y-D)dA= f E « x (y - D) dA. 

Ja J A 

= I E, (-e E ) (y-D)dA +f E 3 <e E ) <y - D) dA. 

-'upper Slower 

CERAMICS CERAMICS 

M = W a E a d 31 (t a + t b )V a . 

Q in • lb f 
M = 1.02x10- 3 Vg-yo^j 

, in • lb f 

ACTUAL MEASUREMENT: M = 0.714 x 10~ J V a VOL ~ f - 



MECHANICS OF PIEZOELECTRICS (CONT.) 


The sensor responds to the applied stress. Assuming that the stress at 
the midthickness is the sensor stress a , , we find that the sensor responds to 

the curvature of the beam |p(- . Again, using the modal expansion, we find that 
the sensor voltage is a function of the curvature of the mass-normalized mode 
shape . This measurement is related to the bending strain of the beam which is a 
generalized displacement. 


SENSOR: PLANT TRANSFER FUNCTION 

THE SENSOR SENSES THE MODAL "DISPLACEMENT" 

°s “ ~J (t s + V T' 
v s - 4 9 3 1 t s< t . + t b»T- 

V s ” ~ 2 ”*$ 931 *s + 

' ' 

“2 

f s = SENSOR CALIBRATION FACTOR = 0.75 



C = FACTOR FOR CURVATURE VARIATION 

/ 3 2 <t> ;(x)\ 

V s = a 2 ?j (t)C^— -L-)= a 2 4 j (t) Cj . 

*THE SENSOR SENSES THE MODAL "DISPLACEMENT". 
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EQUATION OF MOTION FOR PIEZOBEAM 


The partial differential equation of motion for the piezobeam is shown. 

It involves the second spatial derivative of the applied moment . The actua- 
tors are modelled as applying a uniform distributed follower moment over part 
of the length of the beam. The applied moment is modelled mathematically using 
a Heaviside Step function to turn the moment on and another Heaviside Step to 
turn it off spatially. Using the standard modal expansion, the modal equations 
are derived. The coupling of the actuator to the modal equations involves the 
difference in slopes of the mass-normalized mode shapes at the ends of the ac- 
tuator ceramics. 


APPLIED MOMENT M 



m(x) 




E(x) l(x) 


a 2 y \ _ a 2 M(x) 
ax 2 / ax 2 


M(x) = a-|V a [h(x - x-j) - h(x - x 2 )] , 

h ( • ) = HEAVISIDE STEP FUNCTION 
MODAL EQUATIONS: 


y(x, t) = £ ij(t) 0 j(x). 

+ = a^D l V a . 

WHERE Dj = [0j(x 2 )-0j(x 1 )], 
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SYSTEM CONFIGURATION 


The top figure is a block diagram of the control system. The actuators and 
sensors are modelled as non-dynamic real constant matrices. The external dis- 
turbances enter through the shaker. The control approach used in these exper- 
iments is called Positive Position Feedback. This technique uses displacement 
measurements to effect vibration suppression. It can be understood by consid- 
ering the scalar case consisting of two equations, one representing the struc- 
ture or mode £. and one representing a tuned control filter 77. The modal dis- 
placement drives the filter, and the filter coordinate is fed back in turn to 
the structure . The Positive Position terminology can be understood from these 
equations . 


SHAKER 



OUTPUT 


COMPENSATOR 

SYSTEM EQUATIONS FOR SISO: 

STRUCTURE: i + 2f + u 2 ? = gw 2 r) + f (t) 

COMPENSATOR: rj + 2f f w f i) + n = « 2 f ? 


u = MODAL FREQUENCY , f = MODAL DAMPING 

g = GAIN FACTOR , f = EXTERNAL FORCE 

oj f - FILTER FREQUENCY , f f = FILTER DAMPING 

POSITIVE POSITION FEEDBACK (PPF) CONTROL 

• POSITION COORDINATE OF THE STRUCTURE IS POSITIVELY 
FED TO THE FILTER 

• POSITION COORDINATE OF THE FILTER IS POSITIVELY FED 
BACK TO THE STRUCTURE 
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PPF SYSTEM STABILITY 


The system matrix equation for the scalar example shows that the coupling 
of the structure to the compensator occurs in the frequency or stiffness ma- 
trix. This is because displacements are used as measurement quantities. A 
stability analysis of the system equation indicates that stability is main- 
tained if the gain g lies between zero and one . In particular , the point on the 
Nyquist plot which determines stability is the point A. If point A lies to the 
right of the origin, stability is maintained. Point A is the point on the locus 
corresponding to zero frequency. Thus the stability criterion is non - dynamic . 
This is characteristic of Positive Position Feedback and accounts for the im- 
proved robust stability of this method. A root locus for the scalar case shows 
how PPF achieves Vibration Suppression. The filter pole moves toward the imag- 
inary axis while the structural pole moves into the left half plane . Thus , the 
structural pole is stabilized. 



THE CONTROL GAIN FACTOR IS IN THE SYSTEM STIFFNESS 
MATRIX, THUS THE TERM "STIFFNESS CONTROL". 




REALS 
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POSITIVE POSITION FILTER DESIGN 


The Positive Position Feedback compensator is composed of tuned filters 
with transfer function shown below. A simple analog filter realization with 
the desired transfer function is shown. The frequency and damping ratio is se- 
lected based on the results of the control synthesis . 


LAPLACE TRANSFORM OF FILTER 

2 

“f 

T(s) -~2 2 ~ 

s + 2 J f<J f s + cOj 



r f =^ f (R 1 + R 2 )C 2 - 

R 1 = R 2 s 50 K (OHMS) 


THREE MODE CONTROL CIRCUIT 


Several experiments were performed. The first experiments controlled the 
lower modes of the beam individually with one sensor and actuator pair. Then 
the lower modes were controlled together. The control circuit for the control 
of the first three modes is shown along with the component values used. The 
circuit contains more amplifiers than necessary to allow more flexibility in 
the development of the experiment. 



SENSOR 



: | 1 : : 


R6 R7 


± C6 


r® 

-W— I 

f 


R13 

R12 


Rl s 51.5711 
R2= 62.9k 
R3= 70. 2k 
R4= 11 Ok 
R5- 46.4k 
R6= 45.1k 
R7= 51.3k 
R8= 2.92k 
R9= 10.2k 
R10= 2.76k 
R1 1= 994k 
R12= 2.18k 
R13= 10.2k 
R14= 10.0k 


RCTUflTOR 


Cl= !.66uF 
C2= 0.0658uF 
C3= 0.305wF 
C4= 0.0102yF 
C5= 0.137uF 
C6= 0.0059wF 


fil* 741 
A2= INfllOl 


MODE 3 FILTER 


AMPLIFIER 


SISO EXPERIMENTAL RESULTS 


The frequency response functions for the experiments using one actuator 
and sensor pair are shown. These are measured data. The dashed line represents 
the open-loop response, the solid line represents the closed-loop response. 

It can be seen that the control action on the controlled modes greatly reduces 
their response amplitude. In addition, the spillover into uncontrolled modes 
is always stabilizing. This is characteristic of Positive Position Feedback. 
The open-loop and closed-loop free decay for Mode 1 under Mode 1 control is also 
shown . 


MODE ONE CONTROL 



823 







SISO EXPERIMENTAL RESULTS (CONT.) 


Oscilloscope photographs of the free decay of Mode 1 under three mode con- 
trol are shown. Each photograph shows the open-loop decay as the outer enve- 
lope, and the closed-loop decay as the inner trace. Both photographs are of 
the same response only at two different time scales. The settling time of the 
first mode was reduce from about one minute to about one second. 


OPEN AND CLOSED LOOP FREE DECAY OF MODE ONE 
FOR THREE MODE CONTROL 
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SISO EXPERIMENTAL RESULTS (CONT.) 


These tables summarize the open-loop and closed-loop performance for the 
single-input-single-output experiments. Three quantities of interest are 
compared: , and . The first quantity, £, is the damping ratio and is 

a general measure of modal damping. The second quantity, , is inversely re- 
lated to the settling time. The third quantity, , is inversely related to 
the steady-state amplitude of response to sinusoidal excitation ignoring the 
effects of mode shape changes. Depending on the type of dynamic response of 
interest one or another of these quantities is of greater interest. 


EFFECT OF MODE 1 CONTROL ON MODES 1 AND 2 



MODE 1 

MODE 2 1 


? 1 w 1 


f 2 (%) 

h“2 

„ 2 
?2 W 2 

OPEN LOOP 

0.23 

0.0721 

2.27 

0.15 

0.289 

55.5 

CLOSED LOOP 

16.3 

468 

135. 

0.19 

0.366 

70.3 

PERCENT CHANGE* 

7,000 

6.400 

5,800 

26.7 

26.6 

26.7 

PREDICTED^ 

13.3 

3.96 

118. 

0.17 

0.320 

61.7 


EFFECT OF MODE 2 CONTROL ON MODES 1 AND 2 



MODE 1 

MODE 2 ! 



A 

r 2 <*> 


'A 

OPEN LOOP 

0.23 

0.0721 

2.27 

0.15 

0.289 

55.5 

CLOSED LOOP 

0.43 

0.131 

4.00 

12.7 

22.4 

3.95 * 10 3 

PERCENT CHANGE' 

87 

80 

76 

8,400 

7,700 

7.000 

PREDICTED + 

0.37 

0.112 

3.40 

13.3 

23.5 

4.18* 10 3 


EFFECT OF TWO MODE CONTROL ON MODES 1 AND 2 



MODE 1 

MODE 2 | 

?,(%) 

*1 w l 

*1«? 

t 2 (%) 

h^2 


OPEN LOOP 

0 23 

0.0721 

2.27 

0.15 

0.289 

55.5 

CLOSED LOOP 

15.3 

3 81 

94.7 

13.8 

23.7 

4.05 * 10 3 

PERCENT CHANGE’ 

6,600 

5.200 

4.100 

9.100 

8.100 

7 200 

PREDICTED 1 " 

12.7 

3.38 

89.9 

12.7 

228 

4.07 * 10 3 


EFFECT OF THREE MODE CONTROL ON MODES 1, 2, AND 3 



MODE 1 

MODE 2 

MODE 3 | 

?!<%) 

f 1 w 1 

wA 

f 2 (%) 

f 2 w 2 

r J 
f 2 w 2 

f 3 W 

f 3 w 3 

, J 

f 3 w 3 

OPEN LOOP 

0.23 

0.0721 

2.27 

0.15 

0.289 

55.5 

0.27 

1.41 

738. 

CLOSED LOOP 

13.4 

2.63 

51.5 

8.85 

15.7 

2.78 x 10 3 

3.99 

204 

1.04 x 10 4 

PERCENT CHANGE’ 

5.700 

3,500 

2.200 

5,800 

5.300 

4.900 

1,400 

1,300 

1,300 

PREDICTED* 


~ 


" 



~ 

- 

~ 


• PERCENT CHANGE BETWEEN MEASUREO VALUES 
+ PREDICTED CLOSED LOOP VALUES 
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MIMO EXPERIMENT RESULTS 


Two sensors and actuators were used in a multi-input-multi-output exper- 
iment to control the first six modes of the beam. The location of the two sets 
is shown in the figure. The open-loop and closed-loop frequency response func- 
tions are shown for the first eight modes. 


ACTUATOR/SENSOR LOCATIONS FOR MIMO PIEZOBEAM. 
ALL DIMENSIONS IN INCHES. 


1.50 


0.50 


0.25 



1-1.250 -I 


•1 2.50 V—l 



0.0095 

TYP. 


f— 0 020 


OPEN LOOP AND MIMO CLOSED LOOP FREQUENCY 
RESPONSE FUNCTIONS FOR SENSOR 1. 
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MIMO EXPERIMENTAL RESULTS (CONT . ) 


The open-loop and closed-loop free decay responses for the first six modes 
of the beam are shown. 


OPEN LOOP AND CLOSED LOOP FREE DECAY 

MODE 1 MODE 2 
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MIMO EXPERIMENTAL RESULTS (CONT.) 


The two tables summarize the open-loop and closed-loop performance for the 
six mode control case. 

EFFECT OF MIMO CONTROL ON MODES 1 THROUGH 8 



MODE 1 

MODE 2 

MODE 3 

MODE 4 

?,<%> 

fi w i 

. 2 

*1 w 1 

f 2 <%) 

CM 

3 

CM 

Xjr, 

CM CM 

3 

CM 

f 3 <%) 

h w 3 

, 2 
h°>3 

f 4 (%) 

?4“4 

. 2 
?4^ 4 

OPEN LOOP 

0.33 

0.116 

4.11 

0.19 

0.352 

64.1 

0.23 

1.05 

489. 

0.38 

3.54 

3.30 x 10 3 

CLOSED LOOP 

20.0 

5.93 

176. 

24.8 

34.0 

4.65 x 10 3 

8.00 

33.2 

1.38 x 10 4 

4.05 

29.8 

2.19 x 10 4 

PERCENT CHANGE* 

6,000 

5 f 000 

4,200 

13,000 

9,600 

7,200 

3,400 

3,100 

2,700 

970 

740 

560 

PREDICTED^ 

31.5 

8.47 

200. 

18.6 

28.8 

4.44 x 10 3 

13.4 

52.5 

2.06 x 10 4 

5.44 

41.0 

3.08 x 10 4 



MODE 5 

MODE 6 

MODE 7 

MODE 8 

f 5 (% ) 

f 5 w 5 

. 2 
? 5 W 5 

r 6 (%) 

f 6^6 

>• 2 
? 6^6 

r 7 <%> 


„ 2 
f 7 ^ 7 

? 8 (% > 

? 8 W 8 

CM 00 

3 

CO 

OPEN LOOP 

0.39 

5.73 

8.42 x 10 3 

0.37 

7.98 

1.73 x 10 4 

0.34 

10.6 

3.32 x 10 4 

0.36 

11.3 

3.55 x 10 4 

CLOSED LOOP 

0.78 

11.4 

1.68 x 10 4 

0.62 

14.1 

3.14 x 10 4 

0.45 

15.9 

5.09 x 10 4 

0.50 

16.2 

5.20 x 10 4 

PERCENT CHANGE* 

1 00 

100 

100 

70 

80 

80 

46 

50 

53 

40 

43 

46 

PREDICTED^ 

3.24 

46.9 

6.81 x 10 4 

3.03 

64.7 

1.38 x 10 5 

- 

- 

- 

- 

- 

- 


* PERCENT CHANGE BETWEEN MEASURED VALUES + PREDICTED CLOSED LOOP VALUES 
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CONCLUSIONS 


The conclusions drawn from this analysis and these experiments follow. 


• FEASIBILITY OF USING PIEZOELECTRIC MATERIALS AS DUAL-PURPOSE STRUCTURAL 
ELEMENTS/ACTUATORS FOR VIBRATION SUPPRESSION IN LARGE SPACE 
STRUCTURES WAS DEMONSTRATED 

• POSITIVE POSITION FEEDBACK (PPF) AS A VIBRATION SUPPRESSION CONTROL 
STRATEGY WAS IMPLEMENTED 

• USING THE STRAIN SENSOR WHICH MEASURES THE ELASTIC DEFORMATION FOR 
CONTROL WAS SUCCESSFULLY DEMONSTRATED 

• MULTI-MODE VIBRATION SUPPRESSION WAS ACHIEVED WITH DRAMATIC REDUCTION 
IN DYNAMIC RESPONSE 

• NO DESTABILIZING EFFECTS WERE OBSERVED DUE TO EITHER THE SPILLOVER OR THE 
ACTUATOR DYNAMICS 

• A BETTER SYNTHESIS THEORY WHICH PROVIDES PROCEDURES FOR THE SELECTION 
OF GAINS FOR STRONG STRUCTURE/CONTROL COUPLING SHOULD BE DEVELOPED 

• A TRUE ACTIVE MEMBER NEEDS TO BE DEVELOPED AND INCORPORATED INTO MORE 
COMPLICATED EXPERIMENTS 
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